(N 

o 

(N 
+- » 
o 

O 

oo 



From Langevin to generalized Langevin equations 
for the nonequilibrium Rouse model 

Christian Maes 1 and Simi R. Thomas^Rl 
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We consider the Rouse model for the conformational dynamics of an ideal polymer chain subject to steady 
driving. The starting point is a (Markovian) Langevin dynamics of N harmonically coupled particles, some of 
which undergo additional forcing. We compute the effective dynamics for one of the particles by integrating 
out the rest of the chain. The result is a generalized Langevin dynamics for which we give the memory 
and noise kernels, and the effective force. We treat shear-driving as well as a transient regime of constant 
force. These provide natural examples of effective time-evolution starting from nonequilibrium dynamics of 
bath-degrees of freedom. 
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I. INTRODUCTION 

The transfer between different levels of description be- 
longs to the core business of statistical mechanics. That 
is accompanied by the emergence of notions as heat and 
entropy, that are important concepts for bridging the 
gap between microscopic laws and macroscopic behavior. 
And also "new" forces can appear, effectively from inte- 
grating out certain degrees of freedom. These are com- 
monly called entropic or statistical forces and they have 
been discussed since relaxation to equilibrium was first 
considered. Indeed, for equilibrium dynamics, statistical 
forces have their origin in the nature of macroscopic sys- 
tems to evolve towards higher entropy rather than being 
caused by some specific mechanical force. Though much 
systematic progress has been made in the study of these 
entropic forces, there has been less venture into the ef- 
fects of nonequilibrium driving on statistical forces. Only 
over the last decade systematic efforts have been made 
to calculate these nonequilibrium statistical effects. That 
includes Soret-Casimir forces^ but also general consid- 
erations on the law of action and reaction^. 

The first motivation of the present work is to study 
a simple example of statistical forcing emerging from 
integrating out a nonequilibrium environment. Yet, 
the case we treat comes with an extra motivation as 
it opens some questions in the nonequilibrium physics 
of polymers. In contrast to many ongoing studies of 
nonequilibrium polymer rheology, of transport through 
polymers or of mechanical folding and stretching of 
polymers, the present paper considers also steady 
noncquilibria, i.e., where the driving is constant and the 
condition of detailed balance is broken. 
Our working model is the widely studied Rouse model^, 
an ideal chain, where monomers are connected through 
Gaussian springs, and excluded volume effects and 
hydrodynamic interactions are neglected. This model 



holds a special place, as it is the simplest model which 
can be exactly solved to describe phenomena like anoma- 
lous diffusion of polymers in a bath. Moreover, in the 
natural context of polymer melts, which are a collection 
of polymers in a solution, the diffusion of a tagged 
polymer can be described by Rouse dynamics moving 
along a one-dimensional tube embedded in a network 
or mesh of polymers. A question which naturally arises 
in the above context is what the effective dynamics 
of a tagged monomer is when the chain is subjected 
to nonequilibrium driving. We can have in mind that 
the extremal monomers are sheared and that we wish 
to follow a tagged monomer near the middle of the chain. 

Because of the harmonic interaction, the nonequilib- 
rium Rouse model is one of the simplest, yet physically 
interesting examples to understand the effective dynam- 
ics in a driven medium. The equilibrium version of in- 
tegrating out the Rouse model was already carried out 
by D. Panja£. The dynamics of the tagged monomer is 
shown to be non-Markovian with memory kernel having 
a power law decay /i(i) cx i -1 / 2 for short times and ex- 
ponential decay asymptotically in time, 
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The kernel fj,(t) is shown to be the mean relaxation 
response of the polymers to local strain and its behavior 
gives good information on the nature of the diffusion, 
which is anomalous for intermediate times, Ax 2 oc D\[t. 
In this paper we start with the phantom Rouse dynamics 
in the inertial regime and we introduce a nonequilibrium 
driving. The result of integrating out the (other) poly- 
mer degrees of freedom is again a generalized Langevin 
equation (GLE) for the tagged monomer. We show 
that in the overdamped limit, the equilibrium results 
match those of Panja£. We discuss the nonequilibrium 
corrections to the force and memory terms for driving 
of specific nature to obtain some general information 
about statistical forces in nonequilibrium. 
The more systematic and general approach to integrating 
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out degrees of freedom is commonly referred to as the 
Mori-Zwanzig approach^ or the approach via adiabatic 
elimination^ - — . Generalized Langevin equations have 
also been derived in nonstationary environments^ and 
similar in spirit as in the present paper is also the 
generalization where a coarse-graining is added upon 
a coarse-grained description^, or where one Brownian 
particle is described in a nonequilibrium bathi^. In 
the case of nonequilibrium thermostated dynamics a 
generalized Langevin equation has also been derived^. 
We do not follow these general schemes here also because 
we work on the Langevin side of the question, and we 
take no special limits for macroscopic systems or for 
the speed of motion of non-conserved versus conserved 
quantities. Moreover these general approaches are 
less explored for starting with open nonequilibrium 
dynamics for polymers, as they usually rapidly become 
very non-explicit by lack of information about the 
nonequilibrium stationary behavior. 

In the next section we introduce the model and the var- 
ious types of nonequilibrium driving. Sections HVl FVl IVl] 
summarize the method and the results with a discussion 
of the effective dynamical behavior. Finally some more 
essential elements of the computations are collected in 
Section ELD 



II. NONEQUILIBRIUM ROUSE DYNAMICS 

We consider the positions R4, i = 1, . . . , N, of N point 
particles in a three-dimensional domain subject to ge- 
ometric (boundary conditions) and thermodynamic con- 

I 

m f^± = - m f dt'^ N \t - 1>)ki(f) 
at 2 J 



Indeed as an explicit example of a type of Zwanzig 
program—, we will find the validity of a GLE of the form 
(|5|) with rather explicit memory and friction kernels. We 
will discuss the memory kernel n(t) (more generally a 
matrix) , the noise ff(t) and the statistical forcing G (that 
all depend on N) in the cases that we introduce next. 
Obviously, the case of the effective dynamics on another 
coordinate, e.g. the middle one around i = N/2, can 
be reduced to that case. The effective force G can be 
of convolution type and also contain the memory of the 
past trajectory of the tagged monomer. 



straints (like environment temperature and density). The 
particles are harmonically coupled and some are subject 
to further forces, some of which are non-conservative. 
The potential energy is quadratic 

U(R) = ~(Ri-R2) 2 +~{R2-R 3 ) 2 +. ■ .+~{R N -i-Rn) 2 

(2) 

and the force on the ith particle is the sum of systematic 
forces Ki and Langevin forces L^. 

Ki = -ViU + F u Li = -rwyRi + & (3) 

There is an independent standard white noise modeling 
the action of a thermal environment at temperature T 
and friction 7, {£i,a(t)€j,p(tf)) = 2m'ykBT5ij8 a ,pS(t-t'), 
where a and (3 refer to the various directions. The first 
term in Ki is the conservative part of the force. The 
force Fi need not be conservative or constraining and 
will be specified below; that is what we refer to as the 
driving. We then have the equation of motion for the 
time-dependent coordinates Ri(t) 

d R» =• , 

m ~^r = K i + L i ( 4 ) 

with given initial conditions Ri(0), Ri(0) at time t = 0. 
The present paper only considers the linear case where 
the driving Fi is an affinc function of R so that the 
solution can be obtained by diagonalization. We will 
not strictly follow that road and we will not write the 
solution in terms of modes. After integrating out all 
particles but the first one, we obtain explicit information 
about the final equation of the form, 



77177?! (t) + ff") (t) + a (t) + <?W (f ) (5) 

I 

A. Constant forcing 



The simplest case is to assume that the outer end of the 
polymer is being driven under a constant external force 
/. That is a mathematical idealization of a polymer say 
with a charged end, forced under an electric field. As 
there is no confining force for the polymer, that means 
the whole system will move in the direction of the field 
and we discuss that diffusive regime. That is, we take 
U cx t = and Fi = Sjsr t ife x , for some constant field / 
in the x— direction. The simplest example corresponds 
to two linearly coupled degrees of freedom moving in one 
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dimension, with dynamics 

d 2 R d 
m-Tp- = ~k[Ri - R2] - m^—^ + 

m-J^ = -k[R 2 - Ri] - m~fj t R 2 + + f (6) 

for i?i,i?2 £ R. The constant / induces a drift. We 
give here that dimer-case explicitly also because we have 
found that for all finite N (size of original polymer) the 
basic qualitative features of generalized memory and fric- 
tion are unchanged from N = 2, where things are of 
course much simpler. 



B. Shear 



which would replace the dynamics of the x-components 
of the equations (J7J) . 



III. GENERAL METHOD: INDUCTION AND 
RECURRENCE RELATIONS 



Here we imagine the motion of a polymer in a 2- 
dimensional slab of vertical size L in which the outer 
end is subject to shear forcing. 

For explicitness we write out the shearing case again first 
for a polymer of size TV = 2. One monomer is being acted 
on by a shearing force which depends on its y-coordinatc. 
The equation of motion written in Cartesian coordinates 
is then 



md R 2x 

dt 2 
md 2 R2y 

dt 2 
md 2 R lx 

dt 2 
md 2 R\ y 

dt 2 



k[R lx - R 2x ] - m-i 

k[R\y ~ R 2 y] - TO7 



dR 2x 

dt 
dR 2y 

dt 



+ & x {t) + fR 2y 



-k[Ri x - R 2x ] - m 1—^- + £ix(t) 



dR 



.k[R ly -R 2y }- mi ^+^ y (t) 



dt 



(7) 



where / is the shearing speed. 



It is useful here to have a comparison or equilibrium 
reference, where an external potential J7 e xt is added to 
the potential energy U so to trap the outer monomer. In 
other words, again for simplicity of presentation, for the 
case of a dimer, F(t) = —f(R 2 — Q) which derives from 
a confining potential around position Q which holds the 
outer edge of the polymer. 



d 2 R 2 
1 dt 2 
d 2 Ri 



-n{R 2 -R x )- m 7 ^ + £ a (t) - f(R 2 - Q) 
dt 



-K^-R^-mj—l+^t) 
dt 



(8) 



In this section we show the methods and intermediate 
steps involved in reaching our results that will be 
summarized in the next three sections. The general 
method is always to work via iteration and to prove 
results by induction. More precisely, the tagged particle 
equation of motion is directly coupled to a second 
particle which then is coupled to the other N — 2 
particles. If we now assume that first, after integrating 
out these N — 2 particles, the effective dynamics on the 
second particle is of the form ([5]). then we obtain two 
equations: one is the GLE ((3J) (with N there replaced 
by N — 1) and the other is the original equation of 
motion of the tagged particle coupled to the second 
particle. Assuming the structure ((SJ) for TV — 1 with 
specific properties of the memory kernel, noise and force 
constitutes the induction hypothesis. The remaining 
task is then to integrate out that last (second) particle 
and to prove that the induction hypothesis is indeed 
reproduced at size N. The crucial step to discover 
what is the correct induction hypothesis is the case 
N = 2. That is also why the essential first step is to be 
explicit about the case N = 2. We next give more details. 



After integrating out N — 2 particles we arrive at the 
following GLE for the N — 1th monomer, which we label 
with subscript 2 (second particle). (Note that we skip 
vector notation, as we can always reduce the problem to 
more scalar degrees of freedom.) 



d 2 4 N - 
dt 2 



■(*) 



(R 



(iV-l) 



Ri)-i 



dR, 



(N-l) 



dt 



-(*) 



dt'^(t-t')^_(t') 



m m 



(9) 



The tagged monomer R\ is attached to R 2 N 1 ' by a har- monic spring. The force on which is simply given by 

(10) 
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where £1 (t) and £2 {t) are the independent white noise on 
monomers 1 and 2. 

These equations represent a system of 2 monomers, one 
of which is already a coarse grained variable, with mem- 
ory kernel /i^ -1 ' (<), external force G^ -1 ^) and noise 



A second major ingredient in our computation is quite 
naturally to take the Laplace transform of ([9]) and (|10p . 

After integrating out R^ 1 \ we arrive at the following 
GLE for R u 



¥ N \s) = -run / \> +7 2 „ - Ri(0)} - m^sR^s) - R^O)) 

ms^ ly L )(s) + ms-f + nis^ + k 



msjl( N x ) (s) + 771S7 + ms 2 + re 

TOK -(JV-1U \ _l 7 FT - ^2 (Oj-^ll )] 

msfj, 1 -" 1 >{s) + msj + ms + k 



msfl( N - 1 ) (s) + ms7 + ms 2 + re msjl( N 1 ) (s) + TOS7 + ms 2 + re 

&00 (ii) 



That has to be confronted with ([5]) and in particular 
with each of the terms on the right-hand side. In that 
way we obtain recurrence relations for the memory kernel 
jjS N \s), the noise fj^ N \s) and the induced force & N \s) 
on the tagged particle when comparing size N polymers 



with size N — 1: 



re(/x( JV - 1 )( S )+7 + s) 



(msjl( N - 1 ) (s) + 771S7 + ms 2 + re) 



(12) 



GW(«) 



(msjj,( N ^ (s) + TOS7 + ms 2 + re) 



(13) 



rjW(s) = m^W^^^O) - i? l( 0)] + ■ - 2+ ^ f^W 

[rnsfj.^™ L > (s) + 771S7 + s + k) 

(r? (Ar " 1) ( S )+6( S )) (14) 



(ms/i^ -1 ) (s) + TOS7 + ms 2 + re) 



To these we must add "initial" conditions for the recur- 
rence, i.e., to insert the findings for the case N = 2. 
These will enable the correct induction hypothesis. 
Finally, there are the initial conditions to the dynamics; 
the initial conditions (positions and momenta) of all the 
other particles (except the tagged particle) contribute to 
the noise. Their statistical distribution is in principle a 
matter of choice but there are of course dynamically more 
natural choices. We will detail them in Section IVHI 



IV. FREE DIFFUSION UNDER CONSTANT FORCING 

This and the two following sections summarize the 
main results of the logic explained in the previous Sec- 
tion. We always refer to ([5]) for the notation, that we 
have obtained after integrating out all but one of the 



particles. 

A. In general 

Referring to the dynamics ©-(U])-©, we define the 
frequency to as lu 2 = — If us is real (inertial case), 
then the friction kernel //-^ is oscillating with frequency 
u) with decreasing amplitude. When under high friction, 
lu is imaginary (overdamped case), there is monotone 
decay in time. 

In the long time limit we find that the friction kernel 
is always exponentially decaying as fi^ N \t) < e~*/ TjY for 
large times t, where the decay time tn is of the order 
N 2 for large N, see Fig. [T] In the figures are represented 
polymers of different sizes. The discrete data points are 
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Memory vs time (tails) k = 3; \ = 1 



a. 




50 60 70 80 90 100 110 120 



FIG. 1. The long time limit of memory in a semi-log plot. 
The dotted curves are the analytic values and the full curves 
are linear fits. The slope of the various curves is seen to be 
inversely proportional to iV 2 (N is the system size). 



solutions of numerical evaluation of analytical results. 
The lines are fits exhibiting the general nature of these 
solutions, k is the spring constant and \ = m l- Fig. [T] is 
a semi- log plot, which shows the exponential decay of the 
friction kernel in long time. The slopes of various lines 
are proportional to the decay exponent tjv- It is seen in 
Fig. [1] that tjv varies as N~ 2 , which is already known 
to hold in full equilibrium^. We show that this remains 
valid in nonequilibrium. We have also obtained that the 
time of relaxation is in general bounded from below as 
tm > 2/7, see below under Section [VII 11 This is indeed 
clear as the time of relaxation grows with the size of the 
polymer and the tagged particle relaxes fastest when 
it is connected to just one another monomer, tn then 
being 2/7. 

For short times t < r^r, the memory kernel has a power 
law decay; see FigJ5J The figure is a log-log plot of the 
friction kernel versus time. The straight lines exhibit a 
power law decay with exponent —1/2. 

1 

That power law decay of the memory for short times 
agrees well with the anomalous diffusion of polymer o 17 ' 18 
for short times and as also shown by Panja£. Indeed, we 
find in the overdamped limit that 

M W( t ) = 9Mk e -t/Ts (15) 

yt 

where the function C(t) does no longer depend on N 
and remains of order one. Fig[2] shows the behavior for 
different polymer sizes N (for k = 3,x = 1); the curves 
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FIG. 2. The short time behaviour of memory in a log-log plot. 
The dotted curves are the analytic values and the full curves 
are exponential fits. 



are almost overlapping, showing that for short times 
there is no dependence of memory on system size. 

Integrating out the other monomers also creates addi- 
tional coloured noise in the system. The noise rj( N \i) 
is Gaussian with a shifted average and is breaking the 
second fluctuation-dissipation relation transiently. How- 
ever asymptotically, the stationary covariance satisfies 
the second fluctuation-dissipation relation 

vn 

lim (r,(t + r) V (t)) = -/i(r) (16) 

t->oo p 

Finally, the external force / on the outer monomer gives 
rise to a time-dependent statistical force G(t) on the 
tagged monomer and reaches exponentially-fast a lim- 
iting form. The general behavior is exactly similar to 
what we make explicit in the next subsection; see below 
in d20D- 

B. Two monomer case 

To give immediately more explicit formulae we summa- 
rize the results for the dimer-case, which is also used to 
start the recurrence. First the memory kernel (N = 2), 

fi(t) = — e" 7 * /2 [coswt+ ^-sinwt] (17) 

9 

hi '"y 

where uj 2 = — > 

m 4 

For a dimcr, a power law decay for small time is not 
seen, as can be imagined from the fact that a dimer in 
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one dimension has no conformational degrees of freedom. 
The mean squared end-to-end distance (R^ e ) is simply 
equal to the square of the bond length and thus shows 
only simple diffusion. 

The coloured noise is 

ri(t) = -k[Ri[0) - i? 2 (0)][cos(wt) + -^sin(W)]e- 7t/2 

+ k I* sin(u,(t-t')) (t _ t , )/2dt/ 
m J a uj 

It is natural to take the distribution 

p st (R2,R 2 ) = ^ H (19) 

where H = 7 -^- + f (#1 - R 2 f - }R 2 that depends on 
the position i?i of the tagged particle. When averaged 
over that initial distribution (|19[) we get a mean 

(fj(t)) = /[cob(w*) + Bm(wt)]e-rt* 

which is not zero for / ^ 0. Indeed, the monomer R 2 is 
found more on one side than the other due to the force 
/• 

The effective force is found to be 

Git) = f\l - e- 7t/2 (cos(w£) + — sin(wi))] (20) 

exponentially growing to the applied force /. When the 
force would be time dependent, / = f t , the effective force 
gets memory and becomes 

G(t) = k f f s e" 7(t - s)/2 sinw(t - s) ds 
Jo 

That appears to be a general feature of noncquilibrium 
forcing; they create effective forces that themselves de- 
pend on the forcing at all earlier times, we still empha- 
size that the nature of noncquilibrium in this paper and 
in the present example of constant forcing is qualitatively 
different from the case of driving forces discussed in some 
previous worksi^, polymer translocation by an external 
force being one example. Such driving forces are intro- 
duced at the macroscopic level or the level of the GLE 
and hence do not have effect on the nature of the friction 
kernel or of the noise. On the other hand the nonequilib- 
rium in our work is introduced microscopically, such that 
the GLE itself gets modified as a function of the driving. 

C. Limiting cases 

1. The long time limit. Motion after t tn, of 
the tagged particle appears to be diffusive with a 
constant drift /. 



2. Large coupling limit. A large coupling signifies 
that the restoring force between any two monomers 
is very strong and hence the monomers undergo 

high frequency oscillations given by u> ~ 
00. Measurements will typically time-average over 
a few periods. The time-averaged behavior of the 
tagged monomer is again of a Brownian particle 
acted on by a constant force. The effective force 
goes to the constant force /. The time averaged 
total noise goes to white noise and the time 
averaged memory kernel fi(t) disappears. 

3. The overdamped limit. 

The high friction limit refers to the case when the 
viscosity of the medium is so high that the accel- 
eration of the monomers is zero, and the only vari- 
ables which are changing are position. To take the 
overdamped limit in a meaningful manner, together 
with taking the friction coefficient 7 to infinity one 
has to take the mass m of all monomers to zero, 
preserving the product \ = m 7 to be finite. The 
memory kernel then reduces to 

H(t) = ke~ kt/x (21) 

The memory kernel after taking the continuum 
limit is 

= 2^e"^ (22) 

where r = N 2 x/(TT 2 k), as shown before^ under 
equilibrium dynamics. 

V. THE SHEARING CASE 

We are now in two dimensions with shearing at one 
end of the polymer. An interesting result here is that the 
friction kernel fj,( N > (t) for shearing is identical to the case 
of constant forcing (previous section) . To be explicit and 
without loss of essential information we can already state 
the results for N = 2. 

The memory kernel in each direction is given by 

Hix(t) = Vi y (t) = — e~ 7 * /2 [cos(cj£) + -^-sin(wt)] 
m Jo 2w 

which is indeed the same as for a dimer in free space 
under constant forcing. That is due to the fact that the 
external forcing does not couple to velocity but only to 
position. The nature of fi in the two mutually perpen- 
dicular directions and various limits hence remains the 
same. 

The nature of the induced noise however gets modified 
due to the different nature of the external force. The 
noise in general is dependent on the initial positions 
and velocities of all the monomers, and hence picks up 
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additional contributions from the external force which is monomer. Here is the explicit noise function 
coupled to the y component of the position of the first 



Vix(t) = -k[R 1x (0) - R 2x {0)}(cos(ojt) + sin(wi))e~ 7t/2 

id 

+ * [\-^y^ {t - t ' ) \ 2x {t')dt' 

m J to 
+ —k f di'e- 7(t -*' )/2 sm(w(t - t'))[R 2y (0)(cos(ujt') + 

f sinM0) + ^(O)^] 

+ —— [ f dt'dt"e-^ t - t "V 2 sin(u(t - t')) sinfwft' - *"))&«(*") (23) 
tow mw J J 



Putting / = in (|23|) gives us back the noise on a poly- 
mer under constant force (fT5| . The external force couples 
the y component of noise to the dynamics in x direction. 
The initial positions in the y-direction as well as the com- 
ponent of the white noise £, 2y {t) in the y-direction now 
play a role in the dynamics in the x-direction of the sec- 
ond monomer. The y-component of the noise remains 
unaffected by the force, since the external force does not 
couple to the motion in the y-direction. 



Viy(t) = -K[Riy(0) - #23/(0)] (cos ut + 2- sin ^e" 7 */ 2 

Zuj 



of the effective force is 



nR 2 y(0)e 

- fe- 
rn Jo 



_ 7t/2 sin(wi) 

- 7 (i-(0/2^M)) f /. 



)dt' 



(24) 



Now we come to the induced force. The x— component 
I 



G lx (t) = f 



dt' 
Jo 



dt n e - l{t -t")/2 sin ( w ( t _ t ')) 
sm(u{t'-t"))q ly (t") (25) 



Again, as a nonequilibrium effect, the effective force has 
memory. On the other hand its y-component G\ y (t) stays 
zero. Indeed, since the applied force itself is acting in the 
x-direction, there is no reason why the effective dynam- 
ics in the y-direction of the tagged monomer should get 
affected by it. 

Let us now go to the results for general N. The 
memory kernels p, x N \s) and £i y N \s) look as they did in 
the case of constant force (fl"2"|) . with the same "initial" 
conditions. 

The recurrence relation for the coloured noise picks 
up changes due to shearing, similar as discussed un- 
der (f2"B"]) and (f2~4"|); they arise due to the fact that the 
shearing force couples the x and y components of motion. 



fj x N \s) = m^)( S )[4 ff - 1 )(0) - R x (0)} + 



(msfl x N 1 ^ (s) + 771S7 + ms 2 + k) 



(msfl x N (s) + msj + s 2 + k) 



b(N-l) 



R 



(0) 



•7 + Mf _1) W]# 



(0) + tij»-V (0) + ef - 1} (s)/m + W/m} 



[ms 2 + TTtfys + k + msjl x N l '] [ms 2 + mjs + k + msjly" 1> ][ms 2 + mjs + k] 



(26) 



The most interesting aspect of the shearing case is the 
appearance of memory in the induced force. We have 
already seen this in the case of two-monomers (|2"5]) . This 



behaviour persists in general with, 

Fi N \s) = & x N \s)R y N \s) (27) 



8 



ms 2 + 7 s + k + msjlx ^ (s)] [ms 2 + 7 s + k + msjjLy N ^ (s)] 



(28) 



r 



Starting the recurrence with a single monomer where 
g£\s) = /, all subsequent forces can be determined 
using equations (|2"T)) - (|2"5)) . The initial force on a single 

monomer in the y-direction is zero, hence F^ N) (s) = as 
also seen in the two monomer case. 

We studied the asymptotic behaviour of the memory ker- 
nel G x (s) in the same spirit as in lVII II It can be shown 
easily following the same line of arguments that Gx (s) 
decays exponentially in time. In the long time limit, for 
all TV, 

Gf»(t)<e-* 



The memory kernel is given as 

Given the conservative nature of the forces, the second 
fluctuation-dissipation theorem is seen to hold: 



where r = ti — t% and t\ + ti — > 00. 



VI. TRAPPED MONOMER 



VII. DETAILS AND COMPUTATIONS 



Upon introducing an external potential U cx t such that 
the force in the x-direction on the outer edge depends on 
the x-component of the distance of the monomer from a 
fixed origin, the resulting force is not shearing but trap- 
ping. That is thus an equilibrium reference; the force is 
conservative in nature.. The result a rescaling of the fre- 
quency fl 2 = ^f - — \ . The effective force on the tagged 
monomer due to the action of this external potential is 

Fit) = _ J^_ij 1 (t)+-^{i- e -7*/2( C0B n t+ jy si^nt)} 

K + J K + J 211 

(29) 

We recognize the effective spring replacing two springs 
connected in series, 



1 



1 1 



K ef f K f 



This is another way to understand the net restoring force 
on R\ F(t). After all it looks like a trapping potential 
around the origin but of strength K e ff- 
The coloured noise due to unknown initial conditions is 



K 2 



Vi(t) = ^^(i? 2 (0)-i?i(0))e- 7t/2 (cosm+-f-sinm) 
-^-i? 2 (0)e-T t / 2 (cos fit + -J- sin fit) 



n 



m ./n fi 



(30) 



where (vi{t)}*\ = ^e^*/ 2 (cos(fii) + ^ sinfii) for dis- 
tribution 



Let us specify to the case of constant forcing. The 
initial conditions for the recurrence are p,^(s) — 
0; G( 1 )( S )=//s;f/«( S )=0. 

The stationary distribution of ^(O) given qi(0) is 
given by 



Pst = 



-PH. 



where 



ffef-^o)-,,^)] 2 



G 



(JV-l) (JV-l) 

I2 



>(0) 



where we take it that external force has always been on. 
Of course, it is also important here to separate transient 
from stationary behavior. For example, for the force on 
q2, we have at stationarity G^ -1 ' = limt_>oo G^ N ~^ (t) 
and as an illustration we will show using recurrence 
that in the case of constant forcing it always equals the 
originally applied force G^ -1 ^ = /. 

We start from the relation 

lim G {N - 1} (t) = lim sG^-Vts) 

t >oo s >0 

The recurrence relation for the force (fT3")) starts from a 
dimer, 



kGW(s) 



(msft 1 ) (s) + ms7 + ms 2 + n) 



Pat{R%) = 



s(ms7 + ms 2 + k) 
By a simple calculation it is seen that 

lim sF (2) (s) = f 

s — >0 



If now for a polymer of size N — 2 (induction hypothesis) 
lim S F^ N - 2 \s) = f 

s — >0 

We can use the recurrence relation and the property 
lim s >o sp,( N \s) = shown in I VII 11 to see that also 

lim sF^ N - x \s) = f 

s — >0 

as wanted. The stationary distribution thus is given by 

(32) 

The mean noise is 



n , - I e -/8(4C«^" 1) (0)-9l(0) a )-/9S W - 1) (0)) 

Pst - z e 



K 



1. Asymptotic behaviour of memory 

We show here that the memory kernel fl^ N \s) decays 
exponentially in the long time limit, again by recurrence. 



We take the constant force case as simplest example. 
From (jTTj) we see that for a dimcr 



lim e A V 2) (t) = for A < 7 /2 



which translates to 



lims/i (2) (s-A) =0 (33) 



in the Laplace space. 

Let us assume that for a polymer of size N — 1 



lim sfl^ N l '(a — A) = (induction hypothesis) (34) 

s— >0 



and let us choose A = 7/4. That would show that for a 
polymer of size N, 



lim sft,( N \s — 7/4) = lim s 



k(^ n -V ( s _ 7 / 4 ) + 7 + s - 7/4) 



s->0 (m(s - 7/4)//( w - 1 )(s - 7/4) + m(s - 7/4)7 + m(s - 7 /4) 2 + re) 
lim 3s7/4 = 



(35) 



where we have used recurrence relation (|12|) and hypoth- 
esis (f3"4"|) and the fact that limg^o A ( s ) i s a constant, 
as is easy to show. The above result translates to 

lim 

t— ^OO 

Hence in the long time limit, for all TV, 

H (N) (t)<e-^' 2 

which also proves the claim made in IIV Al that the time 
of relaxation is bounded from below by 2/7. 



2. Second fluctuation-dissipation relation in Laplace Space 

The second fluctuation-dissipation relation says that 
the stationary noise auto-correlation function is propor- 
tional to the memory kernel /i(t) through inverse tem- 
perature P, 

(#(t|r)) = ^(r) + 0(i,r) (36) 

such that in the long time limit all terms of the order 
1/t or greater drop out. 

We continue by deriving that relation in Laplace space. 
Let s be the variable in the Laplace space, domain 
|i?e{s}| < 7/2, such that the Laplace transform is well 
defined; see section I VII II 
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pOO poo 

(fi(sW))= / / e-^e-s'*' (r,(tMt'))dtdt' 



Let t' = t + t 



JO 

/oo 
e- s ' t e- s ' T ( ?/ (t)? 7 (t + r))dr 

/OO 
e- s ' T (v{t)v{t + T))dT 
-t 



Let (s + s')i = T (37) 

T dT / e- aV <77(— - -M— - + r))dr 

It- S + S S + S 



a' 



II 



r 

Rewriting (|36p and plugging in the result in (j3"7| 

Mtx^Mtx^ + r)) = ^M(r) + 0(s + ,s') 
s + s' s + s' p 



1 r°° 

(V(s)fj(s')) = —— / e~ T dT 
s t s Jo 



Hence, 

m r 00 _,, 
lim (s + s')(fj(s)fj(s')) = — / e s r ^(r)dr 

= j(£(») + £(-«)) (38) 

is the form of the second fluctuation-dissipation relation 
in Laplace space. 

3. Proof by induction 

We prove our claim that the second fluctuation- 
dissipation relation holds in case of constant force, for 
a polymer of general size N. We give the explicit calcu- 
lation for the case of a dimer. 

Given memory kernel (fTiZj) and noise (fH)) , which is dis- 



/ 777 

e~ ST {-^(T)+0( S + S '))dT 



I 

tributed as p s t(q2, 92) as given in (|32l) . one can calculate 
the correlation function (fj( 2 \s)rj( 2 \s')) s t- Using the re- 
lation 

lim (s + a') = (39) 

s+s'— ¥0 p 

it is shown by a simple calculation that 

lim ( S + S i)(^\ S )^(s')) = ^{fa{s) + &(-*)} 

(40) 

which proves the result. 

To prove it for a general polymer, we use the induction 
hypothesis that for a polymer of size N — 1 the second 
fluctuation-dissipation relation holds: 



lim (a + J)tf»-H')ff S - X) (-0>- - ^{A^W+A^C-)} 

|-s'->0 p 

I 



From the recurrence relations (fT4|i and from (|39|) . one 
easily shows that 
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lim (* + *0<*J (JV) (*)i7 CJV) (*')}* 



K 2 (^ N -^( S )+fL( N - 1 \s)) 



s+s'-^O 



m/3(s/i( JV - 1 )(s) + sj + s 2 + A)(_ s £(iV-i)(_ s ) - s 7 + s 2 + A) 

2 7 k 2 

m/3(s/i( JV - 1 )(.s) + S7 + .S 2 + - ,s 7 + s 2 + A) 
I 



Using the recurrence relations for memory (fT2j). 



K 2 (A (Ar - 1) (s) + /i( w - 1 )(-.5)) 



m 2 (s/2( JV - 1 )(s) + s 7 + s 2 + i)(-s£("-i)(- fl ) - S7 + s 2 + £) 

2k 2 7 

m 2 (s^W( S ) + s 7 + s 2 + A)(_ a £(AD(_ s ) - S7 + s 2 + A) 



which proves the claim. Therefore, the second 
fluctuation-dissipation relation holds for a polymer 
of arbitrary size under the action of a constant force. 

In the shearing case, we are in two dimensions and 
the computations become more involved, but the basic 
recurrence relations remain in place. 



VIII. CONCLUSIONS AND OUTLOOK 

Integrating out degrees of freedom introduces non- 
Markovian noise, effective forces and memory in a tagged 
particle dynamics. That is true in equilibrium as in 
nonequilibrium.. Certain more detailed aspects are also 
unchanged, like the anomalous nature of the memory ker- 
nel for short times which goes into pure diffusion for long 
times, or the TV 2 dependence of the relaxation times. 
In the case of shearing the effective force is no longer 
conservative, and shows important memory effects, and 
the second fluctuation-dissipation relation gets broken. 
These results are not surprising but they are explicit 
thanks to the harmonic interactions, and via iteration 
of recurrence relations. 
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